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Finite-size-scaling ansatz for the helicity modulus of the triangular-lattice 

three-spin interaction model 
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The Berezinskii-Kosterlitz-Thouless-type continuous phase transition observed in the three-spin 
interaction model is discussed. The relevant field theory describes the topological defects involved 
and enables us to perform the renormalization-group analysis. Based on it, we shall propose the 
finite-size-scaling ansatz for the helicity modulus which exhibits the exponent i? — 3/5 for the 
correlation length in the disordered phase. We perform the Monte Carlo simulations to confirm the 
ansatz. Also, we argue its relevance to the ground-state phase transition in the quantum spin chain. 
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I. INTRODUCTION 

The universality observed in the phase transitions is 
one of the most important phenomena to understand the 
interaction effects. For the two-dimensional (2D) criti- 
cal systems other than those with intrinsic anisotropy, it 
is pronouncedly expressed by the conformal symmetry, 
and the correspondingfield theories are characterized by 
the central charge c ^. In the case c < 1, it appears 
to almost specify the universality class, i.e., the possi- 
ble set of the critical exponents [3|- However, for larger 
values of c, there still exist considerable efforts to under- 
stand their universalities, which are possibly related to 
the exotic phase transitions observed in the complicated 
systems. It is widely known that the frustration effects 
sometimes bring about the critical ground states as well 
as the finite-temperature critical points with larger values 
of c > 1. Thus, they have been gathering great attention 
over the years. On another front, the multispin inter- 
actions appear to include the same effect: The exactly 
solved Baxter- Wu model consisting of the three-Ising- 
spin product interaction is the most basic one [3|, which 
belongs to the same universality class as the four-state 
Potts ferromagnet [J . The Ising and the four-state Potts 
criticalities are of c = 1/2 and 1, respectively. Thus, the 
multispin interactions can be expected as another source 
to bring about the larger value of c. 

In this paper, we investigate the three-spin interaction 
model (TSIM) introduced by Alcaraz et al Q. Suppose 
that {k,l,m) denotes three sites at the corners of each 
elementary plaquette of the triangular lattice A (which 
consists of three sublattices Aa, Ab, and Ac), then the 
following reduced Hamiltonian expresses a class of TSIM: 



J 



PH = -77^ Yl ^°^ ('^'= + ^1+ fm) ■ 



kBT 



(1) 



{k,l: 



The angle variables tpk G [0, 27r) are located on sites, 
and the model parameter, the temperature T, will be 
measured in units of J/fce. In the previous paper [6|, 
we discussed an effective field theory for a related model 
(i.e., its clock version) and predicted the one with c = 2, 
which was followed by the numerical confirmation based 



on the transfer-matrix calculations. Here, we perform 
the renormalization-group (RG) analysis to predict the 
Berezinskii-Kosterlitz-Thouless (BKT) 0, [i| type phase 
transition between the critical and the disordered phases. 
In particular, we shall propose the finite-size-scaling 
(FSS) ansatz for the helicity modulus representing the 
stiffness of the corresponding interface model [3,[l3|- The 
discussion goes on in a parallel way with the c = 1 BKT 
transition case [111 ] , and the ansatz includes the exponent 
D for the correlation length. However, since the topologi- 
cal defects involved (for its definition, see below) are not 
described by the scalars (vortexes), but by the vectors 
[3| like the Burgers vectors for the dislocations in the 
2D melting [IJ, [13| , some modifications occur in the RG 
analysis, and then bring about P = 3/5. The FSS ansatz 
permits us to check it directly by the numerical method, 
and then provides solid evidence to support the theoret- 
ical prediction on the instability in the c = 2 criticality. 



II. THEORY 

The fact that the system is invariant under the global 
spin rotations, (pk -^ 'Pk + J2p=ii.h,cJ2i&Api'p^k,i, with 
the condition V'a + i/'b + V'c = (mod 27r), is important. 
Although this U(l)xU(l) symmetry defined by two inde- 
pendent phases, say V'a,b, is not broken and gives the low- 
temperature critical phase, it specifies the possible type 
of perturbations to bring about the high-temperature dis- 
ordered phase. According to Rcf. Q, the following vector 
sine-Gordon Lagrangian density is relevant to the effec- 
tive description: 



c = §-\\dM-) 
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l|N|| = l 



,iN-0(x) 
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The symbol : : denotes the normal ordering and means 
the subtraction of contractions of fields between them; 
II II indicates the norm of the vector, and a gives a short- 
distance cutoff. The parameter j/i stands for the effective 
fugacity to control the appearance of topological defects 
[0|. *(x) is the two component vector field attached at 
the position x in the basal 2D space, so the first term 



represents the interface model, where K gives its stiff- 
ness. The above symmetry is reahzed as the periodicity 
of the field, i.e., $ = $ + 27re„ [a = 1,2), where ei.2 are 
the normalized fundamental vectors of the so-called re- 
peat lattice TZ [IJI isomorphic to the triangular lattice p] . 
Then, the interface can acquire the discontinuity of the 
amount 27rN with N <E 7?.. The second term consists of 
the vertex operators where iKdi^ = eijdj&, and creates 
the shortest discontinuities among possible ones (i.e., the 
length ||N|j = 1). From the RG viewpoint, these topo- 
logical defects controlled by yi arc enough to be kept in 
the theory because these are the most relevant ones (see 
below). 

The reason why we have started with the Lagrangian 
density instead of the vector Coulomb gas (CG) repre- 
sentation Q is that, for the RG analysis we shall em- 
ploy the conformal field theory (CFT) technology, which 
requires the scaling dimensions of local density opera- 
tors and the operator-product-expansion (OPE) coeffi- 
cients among them. While details of their derivations 
on the Gaussian fixed point (i.e., the first term) will be 
given in our future report [l5|, here we shall summarize 
the relevant results to the present case. Since the low- 
temperature critical phase corresponds to the Gaussian 
fixed line parameterized by K, the so-called A4 operator 



x(x)^^lia,*(x) 



(3) 



which shifts the system along the line, is the most impor- 
tant one. Since (A^(x)A^(0))o = (a/?') independently 
of K [r is the distance between and x), it is truly 
marginal. On the other hand, the normalized form of 
the second term in Eq. ^ is given as 
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2/2 = K/2. Thus, W 
E 4) and brings about 



whose dimension is xyj = A'||N 
becomes marginal at K = K* ( 
the transition to the disordered phase for K < K* . For 
later RG argument, the expansion of the operator prod- 
uct >V(x)>V(0) is crucial. While there are 36 terms 
in the double summations with respect to the vector 
charges (say N and N'), the following two cases are 
enough to be taken into account: (i) N -I- N' = and 
(ii) |jN-|-N'|j = 1 (the other terms arc irrelevant here). 
After some calculus using the complex coordinate and 
employing the chiral decomposed form of fields as usual, 
we find that the cases (i) and (ii) mainly give M and W, 
respectively. Then, the expression of the OPE becomes 
as follows 111, [11]: 

>V(x)>V(0) 
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where "• • • " includes the unit operator and the stress 
tensor as well as less singular terms. Consequently, we 
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FIG. 1: (Color online) (a): The contour plot of A^ indi- 
cating the RG flow. Numerical values near the contour lines 
in black denote 10^ x A'^. The straight lines in color show 
j/i(0) = (l-|-f)2/o(0) for some values of f: 0.30 (red dotted line), 
0.04 (pink short-dashed line), —0.03 (green dashed line), and 
—0.20 (blue dot-dashed line), (b) and (c): The scaled cou- 
plings j/o,i/| Ap' ^ are given as functions of the scaled variable 
lAp/^L ' 



can obtain the OPE coefficients CwwM = xw/^/2 and 
Cwww = 2/\/6. The significance of this result is as fol- 
lows: Due to the triangular lattice structure of TZ, three 
vector charges at the angle of 120 degrees to each other 
(e.g., ei, — ei -I- 62, and -62) satisfy the vector charge 
neutrality condition, and then provide the nonzero value 
of Cwww- This brings about the difference from the 
BKT transition (see below). 

Now, the RG equations are derived as follows. Sup- 
pose that the critical fixed point is perturbed by the 
marginal scalar operators C'^(x) (normalized) as £ = 
'^o + ^n^fj-^t^/'^'^^^^' then, the RG equations for the 
change of the cutoff a — > (1 + dl)a are given by 
dXf^/dl = —(1/2) J2i, p C^up^u^p within the one-loop cal- 
culations [131 . For the present case of the couplings j/o 
(= K/K* — 1) and yi, the equations become 



d2;o(0 



-3yl(0^ 



dyi(0 



-2yo{l)vi{l)-Vi{lf. (6) 



dl "^ ' ' dl 

Whereas these are similar to the BKT RG equations [3, 
Q , the y\ term emerges in the second equation due to the 
nonvanishing coefficient Cwww (see also [5|). Now, we 
shall discuss its consequences. Like to the BKT transition 
case, these equations possess a conserved quantity, but 
unlike to the case, it is given as a homogeneous expression 
of degree five in y^ and j/i as 



(yi - yo)' 



y\ + gJ/O 



(AeM). 



(7) 



Figure [Ija) gives the contour plot with arrows to show 
the trajectories of the RG flows. Then, we see that the 
hne yi = yo is the separatrix between the critical and the 
disordered phases, so the small parameter for the phase 
transition can be introduced as j/i(0) = (1 + t)yo{0), with 
\t\ <C 1. Now, with the aid of the conservation law, we 
can analyze the RG flows and obtain the following form: 



yoAl)^\^\^9o,ii\^\^l + rt) 



(8) 



where the superscript "±" indicates the sign of A and 
Tf is a certain function of t given by the initial values of 
the couplings [1^, [201 ■ This form exhibits a kind of self- 
similarity that the trajectories with the initial conditions 
having the same values of t [e.g., all points on the red 
dotted line in Fig. [T]Ja)] fall into a single curve indepen- 
dently of A [the red dotted line (a numerical integration) 
in Fig. [IJb) or [TJ^c)] , and that the t dependence can be 
absorbed by changing the origin of the scaled variable 
according to Tt (the red dotted line overlaps with the 
pink short-dashed line). The explicit forms of gg i are 
complicated while those in the BKT transition are the 
trigonometric or the hyperbolic functions [l^, [20| . Nev- 
ertheless, we can extract some properties because they 
share the basic feature with the BKT transition case. 
On the separatrix t = 0, the solution is simply given by 
yo i(Z) = y{l) = l/3{l + ItJ with It, = l/3y(0). Around 
it,' we find yo,i(0 ^ y{l) - [(1 ± 2)/4] (3/5)^/3 A/y2/3(/) 
when \t\^^^(l/lT,) ^ 1 is satisfied (the upper sign refers 
to the former). Note that these expressions exhibit the 
nonsingularities of yo.i(0 at A = and also that they 
can be regarded as the expansions of the following scal- 
ing forms: 



2/o,i(0 - G + ^T.r '*o,i[A(/ + ?tJ*] 



(9) 



with *o,i(^) — l/i + OiX). To cast these into the forms 
compatible with Eq. ([5]), wc should further replace It, 
by a temperature dependent parameter It = T(/|A|3/^. 
However, since It means the logarithmic scale to obtain 
the initial values yo,i(0) from j/o.i — oo lying almost on 
the separatrix, it is a smooth function and can be rep- 
resented by its value at the transition point, lT,^ near 

Now, we shall estimate the correlation length from Eqs. 
^ and ^. When we write the characteristic logarithmic 
scale as /p satisfying the condition yoilp) = — oo, then 
^(T) ex cxp{If). Figurc^b) exhibits g^ir ~ 1.4) = -cx), 
and A is proportional to T — Tc near T^. Thus, we obtain 
^(T) ex exp[const/(r- Tc)^] with the exponent D = 3/5. 
Unlike to the BKT transition case D = 1/2 0, H], our 
theory predicts the above exponent, which is attributed 
to the nature of the topological defects characterized by 
the vector charges, and more precisely to the nonvanish- 
ing OPE coefficient Cwww. However, our result also 
disagrees with the previous one based on the vector CG 
representation [5|. Therefore, the numerical evidence to 
support our theory is desired. 



In the remainder of the paper, we shall provide the 
evidence by the use of the Monte Carlo (MC) method. 
For this purpose, here we explain the FSS property of the 
helicity modulus to focus on below. The helicity modulus 
is defined as the response of the free energy against the 
long-wave-length (A) twist of the local order field, which 
is proportional to the square of the wave number q = 
27r/A, i.e., T(r) = limg_,od^fiT,q)/dq^ U\M- For the 
present model with the U(l)xU(l) symmetry, the twist 
can be imposed by, for instance, V ^ G Ab; (pi ^ ipi + q-xi 
and V TO € Ac; ipm -^ (pm — Q • x,„, where x; (x,,,) is the 
position vector of the ^th (mth) site and q = q/g is a unit 
vector in the x direction of the basal 2D space. Then, 
one finds the following expression: 



nT{T) 



-H 

2 T 



T\^ sin (^t q • x;^ 



(10) 



{klra) 



where ipt = fk + fi + fm, ^im ^ ^i - ^m, and il is 
the 2D volume of the system. While this is useful for 
numerical estimations of T, its theoretical expression is 
also obtained as follows. According to Ref. 0, the above 
twist corresponds to the following shift in the vector field: 
#(x) — > $(x) -|- q • x(— ei -I- 62), which then brings about 
the increase of the Gaussian part in the free energy den- 
sity {TK/2-K){q^l2). For T < Tc, we can renormalize the 
effects of the topological defects, so that the helicity mod- 
ulus in the full theory is given by T(T) = TK{1 = oo)/2-k. 
As usual, the ratio 7rT(T)/2T exhibits the universal jump 
from 1 to at T = Tc in the thermodynamic limit al- 
though it is rounded off in the finite-size systems. Relat- 
ing to this, Harada and Kawashima performed the MC 
simulations of the 2D quantum XY model, and gave the 
FSS analysis of the helicity modulus to confirm the BKT 
theory being valid there [Tl|. According to their discus- 
sion, here we propose the FSS form for the present phase 
transition: For the finite-size systems with the linear di- 
mension L near Tc, 1+It, and A in Eq. ([9]) are replaced by 
the logarithmic scale ln(T/To) (^o could depend on the 
temperature) and T — Tc, respectively. Then, wc obtain 

^r{T,L) ^_ vl/o{(T-Tc)[ln(T/Lo)]§} ^ ^^^^ 



2r 



ln(L/Lo) 



where Tc and Tq are the parameters to be determined 
based on the FSS ansatz. 



III. NUMERICAL CALCULATIONS 

Here, we shall provide the MC data of the helicity mod- 
ulus and its FSS analysis (the details of our simulations 
will be given elsewhere). Using the standard Metropo- 
lis algorithm, the systems of sizes up to L = 384 were 
treated. The lO'^ MC steps (MCS) were used for sam- 
plings after the 10^ MCS for equilibration. At each tem- 
perature, we performed the independent runs up to 64 
so as to attain reliable statistics of data. In units of 
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FIG. 2: (Color online) Temperature dependence of the he- 
licity modulus. The solid line exhibits 2T/n, and the filled 
circle is T(0) — (,. The inset gives the magnified view near 
the transition temperature, where nT /2T is plotted against 
T. 



fe 



c 




{T-Te)[ln(L/Lo)]-^ 



FIG. 3: (Color online) The FSS plot of the helicity modulus. 
Tc = 0.907 and Lq = hOScaling are employed according to 
the best by-eye scaling. The cross + shows the theoretical 
prediction on the scaling function ^o(O) = 1/3 (see the text). 



the squared lattiee eonstant a^ = 1, the volume is given 
as ri = L'^/C, where C, (= Z/Vs) is the geometric fac- 
tor of A. Then, we obtain the hehcity modulus given 
in Fig. [21 As expected, it exhibits the steep decrease 
and crosses the straight line 2T/tt (see the circles for 
L = 384). Since T stands for the stiffness of the inter- 
face, its reduction to indicates a kind of melting. The 
inset gives the magnified view near the transition tem- 
perature, where ttT/2T is plotted against T. Thus, the 
crossings with the dotted hue provide the finite-size esti- 
mates of Tc, TciL). In Fig. El we show the FSS plot of T. 
According to Eq. pT|) . the parameters Tc and Lq were 
searched until the best collapsing of the scaled data is 



achieved. Actually, we performed it by eye and obtained 
the values, Tc ~ 0.907 and Lq — LOScaling. Then, the 
figure exhibits that the scaling region is narrow in the 
upper side of T^. For L > 192, it continues up to about 
0.25 in the abscissa (see the diamonds and circles). Al- 
though this is due to the limitation on the size treated, 
and further to the property of T showing the jump in 
the thermodynamic limit, we can check the reliability of 
our FSS analysis as follows. As in the case of the BKT 
transition, we can extrapolate the finite-size estimates 
as Tc(L) = rc(oo) + ci/(lnL + ca)^/^, where Tc(oo) and 
ci_2 are the least-squares-fitting parameters. Then, we 
obtain Tc{oo) ~ 0.905, which agrees well with the FSS 
result. Further, as argued, the scaling function is the- 
oretically expected to satisfy ^o(O) = 1/3. We denote 
this condition by the cross in Fig. O and then we find the 
good coincidence with the FSS analysis. With respect to 
the possibility of i/ ^ 3/5, we tried the FSS analysis with 
its value, for instance, 2/5, but the goodness of the scal- 
ing in Fig. 121 could not be reproduced within our search. 
Consequently, these numerical data and the consisten- 
cies among them strongly support our field theoretical 
description on the BKT-type phase transition observed 
inTSIM [H. 



IV. DISCUSSIONS AND SUMMARY 

We shall mention the relevance of our theory to other 
systems. The bilinear-biquadratic spin-1 chain is exactly 
solvable in some cases. The Uimin-Lai-Sutherland (ULS) 
point is one of them, on which the low-energy excitations 
are described by the level-1 SU(3) Wess-Zumino-Witten 
model (c = 2), and by which the massless quadrupole 
and the massive Haldane phases are separated [2^ . This 
ground-state phase transition was clarified to be BKT- 
type, and in terms of the small parameter to control 
the distance from the ULS point, the correlation length 
is given in the identical form as the present one (i.e., 
D = 3/5) [231 . Therefore, we think that the phase transi- 
tion discussed here may share the same fixed-point prop- 
erties with it; we shall discuss this issue more closely in 
our future study [15|. To summarize, based on the vector 
sine-Gordon field theory, we have discussed the continu- 
ous phase transition observed in TSIM; the FSS ansatz 
for the helicity modulus was proposed based on the RG 
analysis. Our theoretical predictions were confirmed by 
the use of the MC simulations. 
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